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The obtainment of the pion structure parameters, in particular its electric and 



^ • magnetic polarizabilities [1], has challenged the experimental and theoretical 

^ community since the early sixties [213] and today remains a field of intense 

activity (see e.g. [1]). Indeed, on the experimental side we expect in the future 
more precise measurements of the pion polarizabilities by the COMPASS col- 
laboration at CERN [5], which probably will help to resolve the discrepancy 
between the dispersion sum rule (DSR) and the chiral perturbation theory 
(xPT) for the difference of the dipole polarizabilities (ai — /?i),r±[SE]- On 
the theoretical side there is also growing interest in the determination of the 
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Fig. 1. Leading- A'^c quark-loop diagrams for the 77 vrvr amplitude. The crossed 
terms are not displayed. 

higher order structure characteristics - the quadrupole polarizabilities of the 
pion. It has been reported that the DSR values of the quadrupole polarizabil- 
ities (^2 ± P2)tv± and (02 + P2)tt0 PISpO] disagree with the present two-loop 
XPT calculations [T1IT2] . 

We contribute to this study by calculating the pion electromagnetic polariz- 
abilities within the Nambu-Jona-Lasinio model (NJL) To our knowledge, 
this is the first dynamical calculation for the quadrupole polarizabilities. The 
NJL model takes into account the quark-antiquark structure of the pion ex- 
plicitly, providing an example of dynamical chiral symmetry breaking in a 
system with nonlinear four-quark interactions. In Ref. [13] the lowest order 
(dipole) polarizabilities have been computed. The alternative approach within 
the large- distance expansion of the extended NJL model has been used in 
[T5] . and for 77 — * tt^tt^ mode in [T6|17fl8j . Here we extend the study of 
Ref. [13] to the quadrupole case, and find a very reasonable agreement with 
the experimental data. For instance, with the empirical values 771^^0 = 136, 
= 93.1 MeV, and with the quark mass M = 300 MeV, we obtain the value 
(0:2 + P2)it0 = —0.144 ■ 10~^ fm^, which agrees within the error bars with the 
value from DSR, {0^2 + ^2)^0 = — 0.171±0.067- 10"'' fm^ [5]. This is an example 
of the general behavior found here: the results of the NJL calculations of the 
dipole and quadrupole polarizabilities are in a good agreement with the DSR 
values. 

Our calculations at the leading- Ai'c order are done according to the Feynman 
diagrams of Fig. [H The analytic expressions for the A and B amplitudes, 
defined below, are derived with the method presented in detail in Ref. [14J. It 
is also necessary to include the l/A/'c-suppressed one-pion-loop contribution, 
as discussed below. The amplitudes are functions of the Mandelstam variables 
related to the 7(^1, ei) +7(^25 ^2) '^"'{Ps) +'n'^{Pi) reaction with the on-shell 
pions and photons. The amplitude 

T{puP2,P3) = e'e^'e^T^., T^. = A{s,t,u)C>^'' + (1) 



is given in terms of the independent Lorentz tensors 
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= - (^^iti^?^'^ + tiP^Ps + niP'zP'i + sp^^pi) , (2) 

where we have disregarded the terms that vanish upon the conditions ei-pi = 
e2-p2 = 0. We have introduced the notation = ^ — ml, ^ = s,t, u. With the 
scalar quantities A and B one obtains the amphtudes = —{A+mlB) and 
H+- = — ml)B for the equal- hehcity and helicity-fiipped photons. The 
general expression for the polarizabilities are conventionally obtained in the 
t-channel. The superscripts i = N,C denote the neutral and charged pions, 
respectively. The dipole, Q.\, (3\, and quadrupole, 0^2,(^2, polarizabilities are by 
definition extracted from the first two coefficients of the Taylor expansion of 
the amphtudes around s = with u = t = ml [19], 



aA\s,t,u) a ( ^ 2 2^ d ■ 2 2^ , \ 
2m. =^X ^ ' ^ ' 7 

= /3i + ^/52 + ---, (3) 

-am.B\s, t, u) = -am. (^B'{0, ml, ml) + s^5'(0, ml, ml) + . . . j 

= («i + A)^ + ^(«2 + /52)^ + ..., (4) 

with a ^ 1/137 denoting the QED fine structure constant. The Born term, 
arising in the case of the charged pion, is removed from the amplitudes. The 
NJL Lagrangian used in this work contains pseudoscalar isovector and scalar 
isoscalar four-quark interactions and is minimally coupled to the electromag- 
netic field. Since we do not take into account explicit vector meson couplings, 
the value of the quark weak coupling constant is = iQ 

The diagrams of Fig. [1] provide polarizabilities which scale as N^. Besides the 
quark one loop diagrams, it is expected that the pion loops will have important 
contributions, mainly in the cases where the tree-level results are absent (in 
the NJL model the chiral counting of meson tree-level results are classified in 
Refs. [Til[20] ). Although formally the pion-loop contributions to polarizabilities 
are at a suppressed level of 1/Nc, the small value of the pion mass entering the 
chiral logarithms enhances them. It is out of the scope of the present work to 
calculate the pion-loop contributions self-consistently within the NJL model, 
which is an involved task even for the basic pion quantities [2T][22] . Instead, 



^ This affects directly the difference of polarizabilities (ai — /3i)'^ which is propor- 
tional to the difference of the low-energy constants Ig — h of xPT, related to the 
TT — > ewy decay. We stress, however, that the expression for {ai — /Ji)*^ evaluated at 
the leading order of the chiral expansion in the NJL model is formally identical to 
the leading order xPT result [1] 
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we shall include the lowest model-independent pion-loop diagram at the 
order, calculated within xPT in Refs. [23|24f25j . and known to be the only 
non-vanishing contribution to the amplitude A at this order in the neutral 
channel. The pion-loop amplitudes are 

= -4^^G.(.), A^, = -2^GAs), (5) 



where 

— 1 °° f s Y 



The pion loop in the charged mode contributes only to the quadrupole polar- 
izabilities, as it starts out with a term linear in s. It has half the strength of 
the neutral quadrupole case. 

The pion loop as well as the a-exchange diagram of Fig. [1] contribute only to 
the amplitude A. Regarding the amplitude B for the neutral (charged) mode, 
only the quark box (quark box + pion exchange) diagrams contribute in that 
case, starting from the order for the dipole polarizabilities and from the 
order for the quadrupole polarizabilities. Thus the combinations {aj + PjY, 
(j = 1,2), to which the B amplitude leads, provide a genuine test of the 
dynamical predictions of the NJL model at leading order in 1/iVc, as they are 
insensitive to the lowest-order xPT corrections. 

All quark one-loop integrals are regularized using the Pauli-Villars prescription 
with one regulator A and two subtractions [26j (in the context of the NJL 
model see also [271128] ). 

[f{M')]pv = f{M') - f{M' + A^) + j^2 df{M' + A') _ 



This procedure is consistent with the requirements of gauge invariance. Here 
M denotes the constituent quark mass, large due to the spontaneous breaking 
of the chiral symmetry. 

It is possible to obtain the results for the quark diagrams in simple analytic 
forms. The expressions are rather lengthy, so we only show a typical result. 
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Table 1 

The NJL model parameters for the charged and neutral channels, with the input 
marked by *. In all cases /* = 93.1 MeV. 



M* fMeVl 


m* \MeV] 


m fMeVl 


G fGeV^l 


A TMeVl 


250 


139 


5.8 


8.49 


964 


250 


136 


5.6 


8.52 


963 


300 


139 


7.5 


13.1 


827 


300 


136 


7.2 


13.1 


827 


350 


139 


8.4 


17.2 


767 


350 


136 


8.1 


17.3 


765 



The factors g-j^ and are the irqq coupling and the pion wave function nor- 
malization, respectively [T¥|l^ . In the case of the box (box + pion exchange) 
diagrams one is dealing with strictly finite integrals, which remain stable upon 
the removal of the regulator, i.e. limA^oo[/(^^)]py = /(^^)- On the contrary, 
for the factors g.,, and the regularization is essential. 

In Table [1] we collect the model parameters, obtained by fitting the physical 
pion mass and the weak decay constant. The parameters of the model are 
the four-quark coupling constant G, the cutoff A, and the current quark mass 
m. These are traded for f^^ = 93.1 MeV, 171^^ = 139 MeV (charged mode) or 
= 136 MeV (neutral mode), and the constituent quark mass, M. 

Our results are displayed in Tables [2l|4l where we show the NJL model pre- 
dictions for three different values of M. In Table [2] we collect also the results 
of the dispersion relations (DR) fit to the Crystal Ball data for the process 
77 — > -K^TT^ [301I3T] (for which the dipole polarizabilities have been taken from 
[32]): the DSR results of Fil'kov and Kashevarov [8], and the xPT predic- 
tions taken from [TT|55] . In Table [3] we include the values of the dipole and 
quadrupole polarizabilities of the charged pion obtained as a result of the DR 
fit [9j to the available experimental data for the total cross section in the 
reation 77 tt+tt" |34ll35|36|37|38|39j in the energy region from the thresh- 
old up to 2.5 GeV. These results are in conformity with the experiments for 
7r~Z ^ 'yn'Z at Serpukhov P0|^ . 7p — 'yir'^n at the Lebedev Phys. Inst. 
[i2] . and at MAMI [43j. The xPT predictions are taken from [12j. 

In Table S] we present the anatomy of our result for the case M = 300 MeV. 
We display separately the several gauge invariant contributions to the polar- 
izabilities: the box (for neutral polarizabilities), box + pion exchange diagram 
(for the charged polarizabilities), the a exchange, and the pion loop ([5]). The 
pion exchange diagram arises only for the charged channel and builds together 
with the box a gauge invariant amplitude. 
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Table 2 

The dipole (in units of 10~^fm'^) and quadrupole (in units of lO^^fm^) neutral pion 
polarizabilities. The first three rows show our NJL model predictions at various 
values of the quark mass, M. 





(ai + (3i)^o 




("2 + fi2)-K0 


("2 - I32)t,0 


M = 250 MeV 
M = 300 MeV 
M = 350 MeV 


1.13 
0.73 
0.50 


-2.05 
-1.56 
-1.32 


-0.33 
-0.14 
-0.07 


46.3 
36.1 
30.6 


DR fit [8] 
DSR [8] 
yPT |33|llj 


0.98 ± 0.03 
0.802 ± 0.035 
1.1 ±0.3 


-1.6 ±2.2 
-3.49 ±2.13 
-1.9 ±0.2 


-0.181 ±0.004 
-0.171 ±0.067 
0.037 ± 0.003 


39.70 ±0.02 
39.72 ±8.01 
37.6 ±3.3 


Table 3 

Same as in Table 


2] for the dipole and quadrupole charged pion polarizabilities. 




(ai + /3i)^± 


(ai - /3i),r± 


("2 ± /32),r± 


(a2 -/32),r± 


M = 250 MeV 
M = 300 MeV 
M = 350 MeV 


0.28 
0.19 
0.13 


10.4 
9.4 
8.3 


0.46 
0.20 
0.10 


23.8 
17.5 
14.2 


DR fit |9] 
DSR [9] 

XPT [33 12] 


0.i»_oQ2 

0.166 ±0.024 
0.16 [0.16] 


13.0l?J 
13.60 ±2.15 
5.7 ± 1.0 [5.5] 


0.133 ±0.015 
0.121 ±0.064 
-0.001 [-0.001] 


25.0l°:« 
25.75 ±7.03 
16.2 [21.6] 



Table 4 

Contribution of various diagrams to the NJL result for the case M = 300 MeV. 





box ± TT-exchange 


cr-exchange 


pion-loop 


total 


(ai ± [3i)^o 


0.73 








0.73 


("1 - /^OttO 


-11.13 


10.57 


-1.0 


-1.56 


(a2 ± P2)-k0 


-0.144 








-0.144 


("2 -(^2)^,0 


5.09 


9.07 


21.97 


36.13 


{ai + [3i)^± 


0.189 








0.189 


(«! - I3i)t,± 


-0.977 


10.36 





9.39 


("2 ± /52),r± 


0.198 








0.198 


(a2 - I32)-K± 


-1.63 


8.87 


10.29 


17.54 



We observe a dependence on the input choice of the constituent quark mass, 
which affects the polarizabihties more severely in the quadrupole than in the 
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dipole case. Let us first consider the sum of tlie polarizabilities, wliicli involves 
only the box contributions. Comparing the sets with M = 250 MeV and 
M = 350 MeV, the box gets reduced in magnitude by about a factor of 2.2 in 
the dipole case and by about a factor of 4.6 in the quadrupole case, both in 
the neutral and charged channels. The faster change in the quadrupole case 
can be understood with a crude estimate: the ratio of the leading-order term 
of the dipole to the quadrupole in the expansion around = is ~ M^. 
This gives the suppression factor 

{a2 + /^2)Af=250 Mev(Q^l + /3i)m=350 MeV ^ ^M=350 MeV ^ 2 (9) 
(a2 + (32)m=350 Mev(ai + /5i)m=250 MeV ^M=250 MeV 



The same happens in the charged mode. The data lies half way between the 
model values for M = 250 MeV and 350 MeV. The best agreement with 
respect to the DR and DSR data is obtained for the set with M = 300 MeV, 
for which we get an overall good description. 

We stress that the values (including the overall signs) of the sums (a2+/32)7r0,7r± 
are totally determined by the gauge-invariant quark box or the box -|- pion 
exchange contribution. This is a straightforward and the most transparent 
result of the presented NJL-model calculation. Moreover, the main part of 
the box contribution comes from the first non-vanishing p^-order term in the 
chiral expansion. Based on this fact we expect that the contact term of the 
3-loop calculation in xPT may also play an important role in reversing the 
signs of the 2-loop order results for these quantities. 

Let us comment on the channels involving the difference of the electric and 
magnetic polarizabilities. i) (oi — f3i)nO'- here the box contribution is largely 
canceled by the scalar exchange. At the p^-order of the chiral counting they 
cancel exactly [T^. The higher-order contributions are quark-mass dependent, 
decreasing as the constituent quark mass increases. The convergence rate is 
slow, at p^-order one reaches only about 50% of the full sum. ii) (ai — i3i)tt±- 
contrary to the neutral channel, the size of the a-exchange diagram for this 
combination is about an order of magnitude larger than the box + pion ex- 
change diagram, and it becomes the most important contribution. The pion 
loops are absent, iii) (02 — ^2)^^'- the pattern observed in ii) repeats itself 
for the quadrupole polarizabilities. However in this case the subleading in the 
1/Nc counting pion-loop diagram has the same magnitude as the cx-exchange 
term. 

In conclusion, we highlight our main results: the NJL model at the lead- 
ing order of the 1/Nc counting yields the right sign and magnitude for the 
quadrupole polarizabilities of the pion, (^2 + /?2)7rO,7r±, being in conformity 
with the DR analysis of the data and the DSR predictions. The sign is sta- 
ble when the model parameters are changed. The magnitude depends on the 
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value chosen for the unobservable constituent quark mass, but the best over- 
all fit to the other empirical data, typically yielding M ~ 300 MeV, also 
yields the optimum values for the polarizabilities. We have also obtained the 
convergence rate for several combinations of the dipole and quadrupole polar- 
izabilities, comparing the chiral order to the full result. We have discussed 
the sensitivity of the different polarizability combinations to the inclusion of 
the one-pion-loop contribution that is non- vanishing at the level of the chi- 
ral counting. The relative contribution of each gauge-invariant combination to 
the polarizabilities has been evaluated and discussed. 
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